Quantum information processing in bosonic lattices 
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We consider a class of models of self-interacting bosons hopping on a lattice. We show that properly tailored 
space-temporal coherent control of the single-body coupling parameters allows for universal quantum computa- 
tion in a given sector of the global Fock space. This general strategy for encoded universality in bosonic systems 
has in principle several candidates for physical implementation. 
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The central problem in quantum information processing 
(QIP) ^ is the ability to control a quantum system in order to 
achieve some predefined purpose like quantum computation 
(QC). In general, a quantum information processor is realized 
by assembling a large number of copies of a given quantum 
system e.g., a qubit, and by making these copies interact in a 
controlled coherent fashion. A crucial issue of any proposal 
for a QIP implementation is its scalability i.e., the realizabil- 
ity, at least in principle, of the above structure for an arbitrary 
large size. An appropriate architecture to achieve this last goal 
is naturally provided by a lattice with sites hosting the pro- 
cessing quantum systems. 

In this paper we shall describe a general scheme for per- 
forming quantum computation with interacting bosonic par- 
ticles in a lattice. The system under study is very general 
and several existing QC proposals could fit into this frame- 
work. These include optical qubits [|2|, [jj Qj, Josephson junc- 
tion qubits [||] and optical lattice loaded with ultracold bosonic 
atoms e.g., Bose-Einstein condensates (BECs) [^, [7|§]. With- 
out focusing on any of these particular implementations, we 
will develop a general framework for encoding qubits and for 
performing universal quantum gates on such encoded qubits. 
Despite their highly practical relevance, we will not discuss 
decoherence issues since they are strongly dependent on the 
specific physical implementation. 

Let us start by casting the problem we are going to ad- 
dress in a more precise control-theoretic fashion. The single 
mode Fock space will be denoted by h := span{|n)}^ . The 
Hamiltonian acting on Ha '•= h® L that we would like to an- 
alyze is given by: 



H(V) 



H.c. 



(1) 



where: (i) A is an index set (the lattice vertices) with L ele- 
ments; (ii) Cj , Cj (j G A) are bosonic creation and annihilation 

operators and ni := c\ci the corresponding occupation num- 
bers; (iii) V := {V}p} x {V[p} C JR x C, is the set of quasi- 
classical 'control' parameters. The Hamiltonian ([j]) represents 
a generalized Bose-Hubbard model [ph, the terms weighted 

(2) 

by the 's account for the non-linear two-body interactions 

whereas the vfp 's are one-body terms describing the hopping 
of the bosonic particles among the lattice sites. The interplay 
of these two terms is known to give rise to a rich quantum 
phase-diagram with insulating and superfluid regions ||(| [Tc| ]. 



The ultimate goal is to find an M-qubit encoding e : 
(<D^)® M i > Ha such that control on the parameters (t) 
in (yj) would enact universal computational capabilities on the 
code. In this paper we will propose such an encoding which 
will enable us to perform universal quantum computation on a 
suitable sector of the Fock space associated with the bosonic 
lattice flTl|], It is worthwhile to keep in mind that even though 
we will consider only the "spatial" interpretation of the single- 
particle modes i.e., spatially localized wave-functions, they 
could even be momentum modes or modes associated to any 
other single-particle wavefunction. 

The qubit. We define the qubit using two lattice sites (dual rail 
encoding). We denote by a,(aj) and bi(b\) the corresponding 
annihilation (creation) operators for the two bosonic modes (i 
is the qubit index). The Hamiltonian of the system has two 
terms: H = Hq + Hi nt . The first term is the sum of all 
single-qubit Hamiltonians: 



H = J2 Hi = J2 £l '* n li + £2 ' 

i i 

+7i,» n a,i + l2,i n b ,i + n{a\ h + di b\) 



(2) 



with n a _i = a- cii, n^i — b\ n is the tunneling rate be- 
tween the a and b modes of the same qubit (/nfra-qubit tunnel- 
ing rate). The second term represents the interaction between 
different qubits: 



Hi n t — 



i¥=3 



a>i a 



a,i ,v a,j 



(3) 



We assume that qubits interact only via the a modes (with /i^ 
the /nfer-qubit tunneling rate; Xij is the Kerr coupling). An 
intuitive picture is given in Fig. [j]. Another possible geometry, 
for example, is to have a common bus to which only the a- 
mode of each qubit is coupled [|T2|l. 

We define the logical state of the qubit by the number of 
particles in the a mode. Thus, n a — defines the logical 
|0)l state and n a = 1 the logical state. The compu- 
tational space is therefore restricted to n a G {0, 1} and we 
will show how to enforce this condition after each gate oper- 
ation. It is important to stress that the subsystems which sup- 
port our qubits are (finite-dimensional subspaces) of bosonic 
modes rather than particles. Therefore, the paradigm we adopt 
here concerning quantum entanglement for systems of indis- 
tinguishable particles is the one advocated in ref. [O]. We 
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FIG. 1 : An example of a qubit array. The a (b) modes correspond to 
the small (big) disks; the b modes are situated in a zigzag geometry 
in order to minimize their interaction. Inter-qubit coupling is given 
only via the a modes. 



describe the main steps of any quantum computation scheme: 
universal set of gates, state preparation and measurement. 

The following set of gates is universal for quantum com- 
putation 0]: {H, P v , Ctt}, where H = ( J ^ 

Hadamard gate, P v = diag(l, e ltp ) is a single-qubit phase 
shift, and C v is a controlled sign flip. We use the more gen- 
eral controlled phase gate C v = diag (1, 1, 1, e ltp ). We enact 
these gates by controlling the time dependence of the param- 
eters characterizing the system. For the single-qubit gates we 
set all fiij = \ij = and we vary only the single-qubit param- 
eters e a ^(t), "f a .i(t) and T,-(t) (a = 1, 2). On the contrary, for 
the two-qubit gate, we keep constant the single-qubit param- 
eters and we control only the inter-qubit tunneling rate /z^ (t) 
or the Kerr coupling Xij (*)• 

Single-qubit gates. In the absence of any external coupling 
([Hj = Xij = 0), the single-qubit Hamiltonian is (for simplic- 
ity we omit the qubit index i): 

Hi =ein 2 a + e 2 n\ + 71 n a + 72 n b + r(a f 6 + a& f ) (4) 

The total particle number n = n a + rib is conserved, since 
[H, n] = 0. Thus, the one-qubit Hilbert space splits into a 
direct sum TL = ©^LgTin, as the Hamiltonian (Q) leaves 
invariant the subspaces Tt n with total particle number n — 
const and dim7i„ = n + 1. In view of this decompo- 
sition of the Hilbert space, we can relabel the Fock states 
as |n a n&) = \n; n a ). Thus, any vector can be written 
as W) — SnSiLo ",' l n ' Since our initial state will 
be a Fock state and since the Hamiltonian conserves the 
total particle number n, the single-qubit wave-function at 
any time will always remain in the subspace TC n , \ip(t)) = 
Sr=o c ™,i(^)l n ) *) G %n- Basically, the only degree of free- 
dom left is the number of particles i (= n a ) in the a-mode. 

For r = 0, the Fock states \n; i) are are also energy eigen- 
vectors, with eigenvalues given by (we omit the label n since 
the total number of particle is conserved): 



Ei 



(n;i\H \n; i) = eii 2 +e 2 (n — i) 2 +j 1 i+y 2 (n-i) (5) 



i = . . . n. We also want our qubit states (defined as before 
|0)i_ = |n; 0), 1 1) i_ = |n; 1)) to correspond to the degenerate 



ground state Eq — E\. This implies the following energy 
degeneracy condition: 



£1 - (2n - l)e 2 + 71 - 72 = 



(6) 



For the Hadamard gate, we keep £fc(t), 7fc(i) constant (and 
satisfying the degeneracy condition (|g)) and we allow only 
a time-dependent tunneling rate r(t). Since the energy gap 
between the ground and the first excited state is AE = 
E 2 - Ei = 4ne 2 + 2(72 - 71) = 2{e\ + e 2 ), we can 
treat the system in a first approximation as a degenerate two- 
level system, ignoring higher level transitions. The time 
evolution (up to a phase) is given by the operator U (t) = 

/' t T(t')dt' 
r , equivalent to a rotation around the a;-axis 

R x ((9) = e~ l6a * . Then we can obtain the Hadamard gate 

as H = R x (tt/4) P x/2 . Similarly, we have NOT = 

iR x (?r/2). In order to confirm this simple analysis, we 

have performed a full time-dependent simulation in the whole 

Hilbert space. We numerically integrate the Hamiltonian (Q) 

for n — 30. If we adiabatically switch the tunneling rate r(t), 

we can control the population of higher levels (and hence the 

leakage form the computational space) to be negligible (in our 

simulation, this is less than for a Gaussian pulse shape 

r(t)). The time scale required for performing single qubit 

gates is about one order of magnitude smaller than the one 

necessary for the two-qubit gate (see Fig. ||). 

To enact the phase shift gate P v we keep the (intra-qubit) 

tunneling zero (therefore we always stay in the computational 

space) and we allow only a time dependence for (some of) 

the other parameters £k(t), 7fe(i). k — 1,2. Since r = 0, 

the time-dependent Hamiltonian is diagonal and we can solve 

the model analytically. Let [0, T] be the time interval during 

which the gate acts. The time evolution of a Fock state is 



(h = 1): \n; i) 



- lTE > |n; i), with E t = ± / Ei(t)dt 



the average value of Ei (t) during the gate operation. Then the 
gate action on the basis states is: 



U v = e 



-iTEa 



diag(l,e-^) 



(7) 



with ip = T(E\ — Eq). In order to have <p 7^ we need 
to violate the energy degeneracy condition (^) by varying any 
of the four parameters £fc(t), 7fc(i). From an experimental 
point of view, the self-interactions Sk might be harder to con- 
trol, since they are related to the collision rates (in a BEC, 
for example). On the other hand, 7^ are related to the energy 
offset of the trapping potential and are conceivably easier to 
control. Thus, we can keep constant any three of these pa- 
rameters and control only the time variation of the remain- 
ing one (say 71). This method gives us considerable free- 
dom in choosing the shape and duration of the pulses 7! (i) 
(the function 71 (t) is not even necessary to be continuous, it 
should be only integrable). Basically, the only condition is 
7 1 (0) = 71 (T) = 72 - £1 + (2n - l)e 2 , such that the two- 
qubit states are again degenerate after the gate; this ensures 
that the phase difference between |0)l and 1 1) |_ is 'frozen'. 

It is important to note that both rotation angles and ip 
characterizing the single qubit gates depend only on the av- 
erage values of r(t) and 7(f), respectively, and therefore they 
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are relatively robust under small fluctuations of the control pa- 
rameters (but they vary linearly with the gate time). 
Two-qubit gate. An important question is: What type of in- 
teractions, together with the one-qubit gates discussed previ- 
ously, are universal? We will discuss two kind of couplings, 
both nonlinear, which achieve this. 



implies the following transformation for the basis states: 

|01> -> (-l)" ll e- l£T |01) (11) 

and similarly for |10). 

(c) 1 11). For n = 2, the Hamiltonian is: 



K 



Xiji 



This is the well-known Kerr Hamil- 



tonian and is used for optical qubits to enact C v [||]. How- 
ever, in usual materials the nonlinearity (the so-called x^) 
is a few orders of magnitudes smaller than what is needed, 
and hence this scheme for producing the two-qubit gate is 
impractical. The Hamiltonian can be easily integrated and 
the gate action on a two-qubit state is simply given by U — 
diag (1, 1, 1, e~ lTxij ), since in our dual rail encoding we al- 
ways have n a ^,n a j £ {0, 1}. We note that, by considering 
excitons in semiconductor quantum dots as bosons jl5[], this 
nonlinearity is the one used to enact the two-qubit gate in the 
QIP proposal of Biolatti et al. [[l(|. 

In the following we analyze in more detail a second univer- 
sal (along with the one-qubit gates) Hamiltonian. 
(ii) Hij = e(nlj+nlj)+7(n a ,i+naj)+^ij(al aj + ai a]). 
This is the Hamiltonian of two-qubits i and j interacting via 
the a-modes. It is identical to the one-qubit Hamiltonian (Q), 
with e\ t i = ei.j — e and 71^ = 7^ = 7, but now we also 



4e y/2mj{t) 



H£>(t) = \ V2m 3 (t) 2s V5/*ii(t) ) (12) 
\/2>y(t) 4e 

In general we cannot integrate this analytically. If fcj (t) 
const, the exact time evolution is 



-3iet 



11) [cOSCi^tH SUlLU 2 t 

UJ 2 

-!^^sinc^(|02) + |20))l (13) 



with 0J2 = \J £ 2 + Again, since we want to recover the 
1 11) state after the gate operation, we impose the condition 
sin(w 2 7 1 ) = 0, hence u 2 T = m 2 Tr, m 2 £ 7L. In this case, the 



evolution of the state is 111) 



-l) r 



-3ieT\ 



11). Together 



have n a> i, n a ,j £ {0, 1}. Since the total number of particles is with the previous condition (i.e., ui{F = m-yir), we obtain 
conserved [Hij, n a .% + n a ,j] — 0, we can neglect the constant 
term proportional to 7 and rewrite the Hamiltonian as (with 
the obvious notation rii = n a> i, rij = n a j) 



Hij(t) = e{n1 + rVj) + faj (t)(a\ a, + a, aj) 




(14) 



(8) 



Given the Hamiltonian (g|), we want to find the control param- 
eter fiij (t) such that the action of the gate is: 



Modulo single-qubit phases Pf , 9 = 7r(mi — ^/m| — 4mf) 
the gate operation on the basis states \n,i rij) is equivalent to 



\m rij) 



\m rij) 



(9) 



C/(T) = diag(l,l,l,e* 1 ) = C 



(15) 



This ensures that the total particle number for each qubit is 
conserved after the gate operation, i.e., there is no leakage 
from the computational space. Of course, during the gate op- 
eration this is not true, since intermediate states like 1 02) do 
not correspond to any logical state, but we will cancel these 
unwanted states dynamically. Again, let [0, T] be the time 
interval during which the gate acts. There are three possible 
cases, depending on the initial state. Since the logical state 
\riin,j)i is the same as the Fock state \niUj), we can omit 
the subscript L (keeping in mind that some intermediate states 
will not correspond to any logical state). The possible input 
states belong to different representations of with total par- 
ticle number n = 0, 1, 2, respectively (n = m + rij). 
(a) 1 00). This case is trivial, |00) -> |00). 



with (fin = 7r(m2 — \frri2 — 4mf). In Fig. |we present a 
full time-dependent simulation for the evolution of the state 
1 11) with mi = 2 and m 2 — 6. We choose a step function for 
the tunneling rate Hij{t). The simulation is in good agreement 
with the exact solution for constant tunneling presented above. 
After extracting the dynamical er %et phase for each qubit, the 
|11) state picks up aphase (fin/n = 6 — 2\/5 1.53, whereas 
the 1 01) and 1 10) states remain phaseless. 

There is an important point to note here: the nonlinear 
term e(nf + nj) in (||) is essential for enacting the gate. If 
e = (or if we replace it with a linear one 7, n, + 7j rij), 
it can be shown that the gate would be equivalent to 11. This 
result is not surprising. The single-body operators entering 



■1 



(b) |01) and |10). For n = 1, we have H^'it) 
Hij(t)cr x (the superscript (1) refers to the total particle num- 



our basic Hamiltonian (|l]) i.e., {c|c 3 '}|'.-_ 1 span, by commu- 

J, 



ber n). Since [H^ (ti), H^ L> (t 2 )} = at all times, we can 
analytically integrate the time evolution to obtain: 

U (1) (T) = e- leT (lcosLUiT -ia x smuJiT) (10) 

with loi = /Ijj = y / Q T Hij (t)dt. Imposing condition (^), we 
require sin(wiT) = 0. Therefore u>iT = miir, mi £ 7L. This 



tation, a Al — u(L) Lie algebra having N := J2i=i c l c i as 
central element. The lattice Fock space splits into (^f^ 1 )- 
dimensional invariant sectors of Al labelled by the eigenval- 
ues N of N i.e., by the total number of bosons. In order to 
universally manipulate M encoded qubits by using just the 
yt^-elements one has to use L ~ 2 M lattice sites. On the 
other hand our encoding scales linearly with the lattice size 
i.e., M = L/2. Thus, the nonlinear term in the Hamiltonian 
(fil) provides an exponential reduction of resources. 
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FIG. 2: Amplitude and phase (in n units) of |11) state during the C, 
gate operation. The dynamical phase e~ 2lEt has been subtracted. 



Preparation and measurement. It is enough to prepare the 
1 . . . 0) l state of the qubit array. We start by preparing two 
linear optical lattices in which the 6-mode bosons are held in 
a zigzag fashion in order to minimize their interaction (see 
Fig. [j]). The next step is to create the middle row in Fig. [l] 
where the a-modes for all qubits will be held. This can be 
done by engineering the confining potential in order to create 
a second minimum for the a-modes. At this stage there is 
no tunneling between any of these wells, Ti = pij = 0, and 
therefore all qubits are in the |0)i_ state (n a ,i = 0). Another 
possibility is to start from a Mott insulator phase, in which 
exact numbers of atoms are localized at individual lattice sites; 
this has been recently demonstrated experimentally Jl0[ ] . 

The measurement technique is conceptually simple - we 
have to detect, for each qubit, the presence or the absence of 
one boson in the a-mode. For an optical lattice, this can be 
done by fluorescence: an atom present will fluoresce under 



the right laser illumination. The middle row of the qubit array 
will be a succession of dark (bright) spots, i.e., the atom is 
absent (present) in the a-mode, corresponding to qubit in state 
|0)l(|1)l). 

In conclusion, we have provided a further example of the 
paradigm of the so-called encoded-universality Q17[], A lim- 
ited i.e., non-universal set of controllable interactions can still 
provide a full computational power in a suitable encoding sub- 
space. We have presented a general framework for performing 
encoded universal QC on systems of self-interacting bosonic 
particles hopping on a lattice. Our strategy requires the abil- 
ity to control in space and time the one-body couplings of the 
system. A summary of the parameter dependence for the gate 
operations is shown in Table |l] (only one of the two nonlinear 
interactions faj and Xij are sufficient and therefore they can 
be used alternatively, depending on the system). Possible im- 
plementations of this scheme include optical qubits, Joseph- 
son junctions and BEC in optical lattices. 
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TABLE I: A minimal example of time dependence for the control 
parameters. The last line is the degeneracy condition fl6J). 
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